function [ SX, mu, st ] = featureScale( X )

[M, N]  = size(X); 
mu = mean(X);
st = std(X);

SX = (X - repmat(mu,M,1)) ./ repmat(st,M,1);

end

